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The primordial inflation dilutes all matter except the quantum fluctuations which we see in the 
cosmic microwave background (CMB) radiation. Therefore the last phases of inflation must be 
embedded within a beyond the Standard Model (SM) sector where the inflaton can directly excite 
the SM quarks and leptons. In this paper we consider two inflaton candidates LLe and udd whose 
decay can naturally excite all the relevant degrees of freedom besides thermalizing the lightest 
supersymmetric particle (LSP) during and after reheating. In particular, we present the regions of 
the parameter space which can yield successful inflation with the right temperature anisotropy in 
the CMB, the observed relic density for the neutralino LSP, and the recent Higgs mass constraints 
from LHC within the MSSM with non-universal Higgs masses - referred to as the NUHM2 model. 
We found that in most scenarios, the LSP seems strongly mass degenerated with the next to lightest 
LSP (NLSP) and the branching ratio B s — >> very close to the present bound, thus leading 

to falsifiable predictions. Also the dark matter interactions with XENON nuclei would fall within 
the projected range for the XENON1T experiment. In the case of a positive signal of low scale 
supersymmetry at the LHC, one would be able to potentially pin down the inflaton mass by using 
the associated values for the mass of the stau, the stop and the neutralino. 



I. INTRODUCTION 

The primordial inflation must explain the seed per- 
turbations for the cosmic microwave background (CMB) 
radiation [1 j, and after the end of inflation the coher- 
ent oscillations of the inflation must excite the Standard 
Model (SM) quarks and leptons at temperatures suffi- 
ciently high to realize baryons and dark matter in the 
current universe [2j[3]. In this respect, it is vital that the 
last phase of primordial inflation must end in a vacuum 
of Beyond the Standard Model (BSM) physics which can 
solely excite the relevant degrees of freedom required for 
the success of Big Bang Nucleosynthesis (BBN) , see for 
a review [4]. 

Inflation needs a potential which remains sufficiently 
flat along which the slow-roll inflation can take place in 
order to generate the observed temperature anisotropy in 
the CMB. The low scale supersymmetry (SUSY) guaran- 
tees the flatness of such flat directions at a perturbative 
and a non-perturbative level (for a review see [5 ), besides 
providing a falsifiable framework for the BSM physics, 
see [6]. Furthermore, the lightest SUSY particle can be 
absolutely stable under R-parity, and thus provides an 
ideal cold dark matter candidate [7 . 

The flat directions of SUSY, especially the MSSM pro- 
vides nearly 300 gauge- invariant F-and D-flat directions 
[H[9], which are all charged under the SM gauge group. 
Out of these flat directions, there are particularly 2 D- 
flat directions: udd and LLe, which carry the SM charges 
and can be ideal inflaton candidates [T0HT2] . Here 2, d 
correspond to the right handed squarks, L corresponds 
to the left handed slepton, and e corresponds to the 
right handed (charged) leptons. Both the inflaton can- 



didates provide an inflection point in their respective 
potentials where inflation can be driven for sufficiently 
large e- foldings of inflation to explain the current universe 
and explain the seed perturbations for the temperature 
anisotropy in the CMB [TU1 fT2] . 

The inflaton in this case only decays into the MSSM 
degrees of freedom which thermalize the universe with 
a temperature, Tr ~ 10 8 GeV [T3] . This temperature 
is sufficient to excite the degrees of freedom which are 
needed for the lightest supersymmetric particle (LSP) to 
get a relic density that matches observations. It is then 
natural to ask whether there exists any parameter space, 
where both successful inflation and thermal dark matter 
abundance can be explained simultaneously [I4j [15] 1 . 

Recently both the ATLAS and CMS experiments have 
recorded hints of a Higgs boson. With an integrated lu- 
minosity of respectively 4.8 and 4.9 fb _1 and a centre of 
mass energy of ^/s = 7 TeV the ATLAS experiment re- 
ported an excess of events in the H — >> ZZ* — >• 4/ (where 
I is either electrons or muons) and the H — >• 77 channels 
[18j[l9]. In addition, ATLAS published a broad excess 
in the H — >• WW* — » Ivl'v channel from a combined 



1 Inflationary models embedded within a hidden sector with SM 
gauge singlets suffer a serious drawback - it is not at all clear why 
and how such an inflaton would decay solely into the SM degrees 
of freedom. A hidden sector inflaton can couple to many other 
hidden sectors of the BSM, given the fact that any stringy con- 
struction of BSM produces large number of hidden sectors within 
landscape, which can in principle accommodate the Kaluza-Klein 
dark matter as a candidate |16| . The top-down construction of 
inflation generically excites the hidden sectors predominantly as 
compared to the visible sector fields, see \17\ . 
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analysis with these two channels [63]. All these signals 
would point out towards a ~ 125 GeV Higgs boson while 
also excluding masses outside the [116,131] GeV range (a 
part from a small window between 238 and 251 GeV). In 
addition to these results, CMS also observed an excess of 
events but pointing towards a Higgs mass of ~ 119 GeV, 
which was not observed by the ATLAS collaboration. 

In this paper we explore scenarios which lead to a Higgs 
mass in the allowed range, a correct relic density for 
the neutralino LSP and inflaton properties which are in 
agreement with the CMB data. We use a variant of the 
MSSM, which differs from the constrained MSSM and 
minimal supergravity model (mSUGRA) [22]. Recently, 
it was shown indeed that the specific value of rrth — 125 
GeV is hard to accommodate within mSUGRA which 
has a rather restrictive parameter space with degener- 
ate scalar masses mo, and gaugino masses m^, at the 
grand unification (GUT) scale [23-25 . For high values 
of mo and it is possible to get mSUGRA regions 

with right dark matter abundance [26j [27] , but they may 
not be accessible at the LHC scale. It was pointed out 
nevertheless that a variant of MSSM with non-universal 
Higgs masses [28] known as NUHM2 can accommodate 
both the 119 and 125 GeV Higgs mass values and the 
observed relic density for the LSP. It is then natural to 
ask whether, in this framework, the same set of param- 
eters leads to an inflaton mass which is compatible with 
the CMB observations and the current particle physics 
spectrum at low energies. 

Since in our case the inflaton candidates are gauge in- 
variant, by using the renormalization group equations 
(RGE) at one loop level, one can evaluate the mass of 
the inflaton, m^, from the scale of inflation to the scale 
of LHC. This eventually will enable us to relate the infla- 
ton mass with the dark matter parameter space and the 
CP-even Higgs mass. 

The plan of this paper is to briefly discuss in Section 
II the properties of inflation and various observables of 
CMB, which can be satisfied by the two flat direction 
candidates, udd and LLe. We will also discuss how these 
two candidates can also generate departure from random 
gaussian fluctuations which can be verified by the forth- 
coming satellite experiment Planck [31]. In section III, 
we identify benchmark points where the constraints on 
inflation and the LSP relic density converge towards a 
Higgs mass of about 119 GeV or 125 GeV. In section 
IV, we perform a broader scan of the NUHM2 param- 
eters in order to delineate the regions of the parameter 
space where the neutralino has the observed relic den- 
sity and the Higgs mass falls within the allowed range. 
In the final section we will discuss how LHC observables 
and also direct dark matter searches (in particular the 
XENON1T [58] experiment) can be used to probe the 
NUHM2 parameter space and eventually pin down the 
inflaton mass. 
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FIG. 1: A schematic drawing of inflationary potential for ei- 
ther udd or LLe as shown in Eq. ^ . Inflation happens near 
the inflection point as shown by ^inflation = </>o , inflation ends 
at the point of enhanced gauge symmetry, where the entire 
(MS)SM gauge symmetry is recovered. The physical mass 
and couplings at hi gh sca le 0o and 4>lhc are related via RGE 
described by Eqs. fll9|20 ). 



II. INFLATION, CMB OBSERVABLES & 
RENORMALIZATION GROUP EQUATIONS 



A. Inflaton Candidates: Flat Directions of Squarks 
and Sleptons 



In Refs. [TOl [TTJ [I4j [15] the authors have recognized 
two D-flat directions which^can be the ideal inflaton can- 
didates, because both udd and LLe flat directions are 
lifted by higher order superpotential terms of the follow- 
ing form which would provide non-vanishing A-term in 
the potential even at large VEVs 2 : 



(1) 



2 Note that the R-parity is still conserved, both the superpo- 
tentials: udd and LLe do not appear at the renormalizable 
level, they are instead lifted by non-renormalizable operators, 
see Refs. [10] 1111 1141 fT5] for a detailed discussion. Further note 
that both the operators vanish in the vacuum which is shown as 
4> = in Fig. [I] Rest of the MSSM flat directions are lifted by 
hybrid operators of type W D $ n_1 ^/Mp -3 where <E> and ^ are 
MSSM superfields, see 5 8 9 . Such operators do not induce 
non-renormalizable A-term during inflation which is relevant for 
inflation, see Refs. 10 11 for discussion. 
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where A ~ 0(1) 3 . The scalar component of superfield, 
denoted by (/>, is given by 4 



u + d + d 



L + L + e 
y/3 ' 



(2) 



for the Sdd and LLe flat directions respectively. After 
minimizing the potential along the angular direction 
(<3> = $e lB \ we can situate the real part of (j) by rotating 
it to the corresponding angles # m in- The scalar potential 
is then found to be [TU1 ITT] 



(3) 



where and A are the soft breaking mass and the A- 
term respectively {A is a positive quantity since its phase 
is absorbed by a redefinition of during the process) 5 . 
The masses for LLe and udd are given by: 



9 
L 



- m~ 



u a d 



(4) 
(5) 



These masses are now VEV dependent, i.e. m 2 ((/>). The 
inflationary perturbations will be able to constrain the 
inflaton mass only at the scale of inflation, i.e. 0o, while 
LHC will be able to constrain the masses at the LHC 
scale. However both the physical quantities are related 
to each other via RGE as we will discuss below. For 



where a 2 <C 1, there exists a point of inflection (0 O ) in 
V{<t>) 6 , where 



AvTo 



0{a 2 )., 



f"(0o) = O, 



at which 



V(<j> ) = -m^ + O(a'), 
V'(ct> ) = 4a 2 m 2 ^ + 0(a 4 ) . 

2 

777 

V"'(<Po) = 32-^ + 0{a 2 ) . 



(7) 

(8) 

(9) 
(10) 

(11) 



Prom now on we only keep the leading order terms in all 
expressions. Note that inflation occurs within an inter- 
val 7 



|0-0O 



60M|' 



(12) 



in the vicinity of the point of inflection, within which 
the slow roll parameters e = (Mp/2)(V /V) 2 and 77 = 
Mp(V ff /V) are smaller than 1. The Hubble expansion 
rate during inflation is given by 

1 777000 



inf 



45 Mp 



(13) 



In order to obtain the flat potential, it is crucial that 
the A((/> )-term ought to be close to ra^^o) in the above 
potential Eq. ([3]). This can be obtained within two par- 
ticular scenarios: 



40777 j 



• Gravity Mediation: 
(6) In gravity-mediated SUSY breaking, the A-term 

and the soft SUSY breaking mass are of the same 
order of magnitude as the gravitino mass, i.e. ~ 
A [35]. 



The exact value of A is irrelevant for the CMB analysis, as it 
does not modify the CMB predictions. However it is possible to 
extract its value by integrating out the heavy degrees of freedom. 
If the origin of these operators arise from either SU (5) or £0(10), 
then the typical value is of order A ~ O(10~ 2 ) for £0(10) and 
A ~ 0(1) for 577(5), as shown in Ref. [14] . 

The representations for the flat directions are given by: uf = 

71^' = ^Tf^ d k = 73^- Here 1 < a '^7 < 3 are color 
indices, and 1 < k < 3 denote the quark families. The^flat- 
ness constraints require that a / /3 / 7 and j / k. = 

° ^ £5 = ^= ( t) , e fc = ^0, where 1 < a, b < 2 



1 

V3 



^3 V3 \ 

are the weak isospin indices and 1 < k < 3 denote the lep- 
ton families. The flatness constraints require that a 7^ b and 
i / j / fe. Note that the cosmological perturbations do not care 
which combination arises, as gravity couples universally. 
Note that the supergravity corrections do not spoil the inflation- 
ary potential. Typically supergravity corrections lead to Hubble 
induced mass corrections, but in our case the Hubble para mete r 
during inflation is always much smaller than m^, see Eq. ( |13| ), 
for a detailed discussion see Refs. 10 Hi. 



• Split SUSY: 
Normally in Split SUSY scenario where the scale 
of SUSY is high and sfermions are very heavy, the 
A-teim is typically protected by R-symmetry, see 



6 The value of a during inflation could be small, i.e. a ~ 10 10 , 
but it runs dynamically from the GUT scale where A 2 = 40m^ 
to the required value at scale of inflation via the RG-equations. 
For a detailed discussion see Ref. |15| . 

7 For a low scale inflation, setting the initial condition is always 
challenging. However in the case of a MSSM or string theory 
landscape where there are many false vacua at high and high 
scales, then it is conceivable that earlier phases of inflation could 
have occurred in those false vacua. This large vacuum energy 
could lift the flat direction condensate either via quantum fluctu- 
ations |32] . however see also the challenges posed by the quantum 
fluctuations |33| . or via classical initial condition which happens 
at the level of background without any problem, see |34| . 
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Refs. [36l [37j , as a result the A- term could be very 
small compared to the soft masses. However, if the 
Yukawa hierarchy arises from the Froggatt-Nielsen 
mechanism, then the A-term can be as large as that 
of the soft mass, i.e. ~ A, as in the case of 
Ref. [38]. 

Keeping low scale and high scale SUSY breaking sce- 
narios in mind here we will consider a large range of 
(m^, 0o ) to match the cosmological observations. 

B. Cosmological Observables 

1. Gaussian fluctuations and tensor to scalar ratio 

The above potential Eq. ([3| has been studied exten- 
sively in Refs. [TTJ [39j [40]. The amplitude of density 
perturbations Sh and the scalar spectral index n s are 
given by: 

_ 8 msMp 1 . 2 a— r- 

Sh = -7= T2^^2 sm [^COBE v A 2 ] , (14) 

V57T <Pq L\ 

and 

n s = l- 4^ cot[AT C oBE\/A2], (15) 
respectively, where 

A 2 ^ 900a 2 ^ BE (^) 4 . (16) 

In the above, A/cobe is the number of e- foldings be- 
tween the time when the observationally relevant per- 
turbations are generated till the end of inflation and fol- 
lows: A/cobe ^ 66.9 + (l/4)ln(F(0 o )/Mp) - 50. Since 
the perturbations are due to a single field, one does not 
expect large non-Gaussianity from this model (/jvx < 1, 
see Ref. [41]). 

In Fig. [2] we have explored a wide range of the infla- 
ton mass, m^, where inflation can explain the observed 
temperature anisotropy in the CMB with the right am- 
plitude, Sh = 1.91 x 10 -5 , and the tilt in the power 
spectrum, 0.934 < n s < 0.988 [1 . Fig (2] represents the 
inflation energy scale versus the mass of the inflaton. The 
configurations which fit the observed values of Sh and n s 
are shown in blue. Although we have restricted ourselves 
to VEV values below the GUT scale, the model does pro- 
vide negligible running in the tilt which is well within the 
observed limit. 

Here we have allowed for a wide range of and </>o 
values because ultimately we want to show that infla- 
tion can happen within low-scale SUSY scenarios from 
high-scale SUSY breaking soft-masses (cf the split-SUSY 
scenario [38] ) . 

In this paper we will mostly consider scenarios where 
the scale of inflation is low enough that one would not 
expect any observed tensor perturbations in any future 




0.5 1 1.5 2 2.5 3 3.5 4 



m^[10 5 GeV] 

FIG. 2: (0o, rucf,) plane in which inflation is in agreement with 
the cosmological observations of the temperature anisotropy 
of the CMB fluctuations. The blue region shows the inflaton 
energy scale and inflaton mass which are compatible with the 
central value of the amplitude of the seed perturbations, Sh = 
1.91 x 10~ 5 , and the 2a allowed range of spectral tilt 0.934 < 
n s < 0.988 1 . Note that we restricted ourselves to inflaton 
VEVs 0o below the GUT scale. 



CMB experiments. To obtain large observable tensor to 
scalar ratio, r, one would have to embed these inflaton 
candidates within N = 1 supergravity (SUGRA). This 
would modify the potential with a large vacuum energy 
density besides providing SUGRA corrections to mass 
and A-term [42j[43]. One could then obtain r ~ 0.05 for 
both inflaton flat directions: udd and LLe as shown in 



2. Non-Gaussianity: interplay between udd and LLe 

Note that within MSSM flat directions, udd and LLe 
are two independent directions which are lifted by them- 
selves. In principle both the flat directions can be lifted 
by higher order terms. One can imagine that either of 
udd or LLe are lifted by higher order superpotential op- 
erators in Eq. ([!]), such as <I> 9 , <I> 12 ,---, see [44], while 
the other is lifted at <I> 6 . This will create a hierarchy in 
potential energies between the two flat directions. One 
will have a large vacuum energy density compared to the 
other. 

This second flat direction which is lifted by lower or- 
der operators is sometimes known as a curvaton in the 
literature, see for a review [2]. Therefore, either udd or 
LLe could be the inflaton or the curvaton depending on 
term in the superpotential which lifts them. In this case 
the curvaton typically slow rolls and decays later, which 
is responsible for generating a sizable non-Gaussianity 
in the squeezed-limit, which is given by [31] and can be 
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detectable by the Planck satellite: 



rsquezed 
JNL 



G(l)h 2/S ~ 0(b) - (9(1000) . (17) 



for h ~ 10 _1 — 1CT 5 , where h denotes the SM Yukawa 
couplings. The /nl depends on a particular decay chan- 
nel of the squarks and sleptons, therefore the Yukawa 
couplings appear in the analysis. In fact the smallest 
Yukawa coupling dominates /jvl- In this scenario the 
perturbations are mainly created by the flat direction 
which is lifted at the lowest order in the superpotential. 
When there are two flat directions which are lifted simul- 
taneously, the thermalization process is delayed a lot due 
to a phenomenon known as kinematical blocking of the 
decay channels [45]. This explains the non-Gaussianity. 
However, as we shall see the rest of our analysis will be in- 
dependent of whether significant non-Gaussianity is gen- 
erated or not after inflation. 



CMB fluctuations probe the inflaton potential at the in- 
flationary scale. At low energies the inflaton properties 
can be probed by the LHC from the masses of the squarks 
and sleptons. 

The inflaton mass and the non-renormalizable A term 
in the inflationary potential are both scale dependent 
quantities, and they can be tracked down to lower en- 
ergies by using the RGE. In [TTJ [I4j [15], it was shown 
that at one loop level for the relevant flat direction, udd: 



dfi 
dA 
dfi 



1 

1 



8 . 



^-(yMs^ + ^M^). 



(19) 



where ft = £lq 
For LLe: 



is the VEV at which inflation occurs. 



3. Reheating, thermalization, & thermal history 

Instant reheating and thermalization occurs when a 
single flat direction is responsible for inflation and struc- 
ture formation. This is due to the gauge couplings of the 
inflaton to gauge/gaugino fields. Within 10 — 20 infla- 
ton oscillations radiation-dominated universe prevails, as 
shown in Ref. [13 . The resultant reheat temperature at 
which all the MSSM degrees of freedom are in thermal 
equilibrium (kinetic and chemical equilibrium) is given 

by m 



T, 



rh 



2 x 10 s GeV. 



(18) 



Since the temperature of the universe is so high, it imme- 
diately thermalizes the LSP provided it has gauge inter- 
actions. The LSP relic density is then given by the Stan- 
dard (thermal) Freeze-out mechanism. In particular, if 
the neutralino is the LSP, its relic density is determined 
by its annihilation and coannihilation rates. 

The advantage of realizing inflation in the visible sec- 
tor is that it is possible to nail down the thermal his- 
tory of the universe precisely. At temperatures below 
10 — 100 GeV there will be no extra degrees of free- 
dom in the thermal bath except that of the SM, therefore 
BBN can proceed without any trouble within low scale 
SUSY scenario. This reheat temperature is marginally 
compatible with the BBN bound for the gravitino mass 
^3/2 > O(TeV). It is also sufficiently high that various 
mechanisms of baryogenesis may be invoked to generate 
the observed baryon asymmetry of the universe. 



/i 



dfi 
dA 
dfi 



1 ,3„ 2 



10 



Mhl) 



(20) 



where Mi, M 2 , M 3 are [7(1), 517(2) and 577(3) gaugino 
masses, which all equate to at the unification scale, 
and #i, #2 and gs are the associated couplings. To solve 
these equations, one needs to take into account of the 
running of the gaugino masses and coupling constants 
which are given by, see [35] : 



= 0, 



(21) 



with a\ = 11/167T 2 , ct2 = 1/167T 2 and ct\ = — 3/167T 2 . So 
every point in (mo,mi/ 2 ) (where mo and mi/ 2 denote the 
scalar masses and the gauginos at the unification scale 
respectively) plane can now be mapped onto (</>o? m (/>) 
plane 8 . 



III. NUHM2 SCENARIO & CONSTRAINING 
THE NUHM2 PARAMETER SPACE 

The NUHM2 is a variant of MSSM with non-universal 
soft breaking masses mi and m2 which are independent 
for both Higgs doublets [5TH53] . The universality of 
scalar masses mo at the unification scale, i.e. GUT scale, 
is still assumed, but in NUHM2 model, they are different 
from mi and m2- It is well-known that the Higgs masses 



C. Renormalization Group Equations 

Since the inflaton carries the SM charges and they are 
fully embedded within the MSSM, it is possible to probe 
various regions of the parameter space for inflation. The 



The RGE equations also exhibit explicitly that the fine-tuning 
required to match and the A-term at the inflection point 
can be obtained from the running of the gauge couplings, see 
Ref. 15 . At the LHC scale the ratio of soft masses and the 
A-term is order one. 
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can be written as, see [52], [53]: 

m\ (1 + tan 2 f3) = M\ tan 2 f3 - /i 2 (tan 2 /3 + 1 - A^ 2) ) 

- (c + 2c M ) tan 2 /3 - A A tan 2 /3 - -m|(l - tan 2 /3) - 

2 (22) 

and 

m|(l + tan 2 /3) = - /x 2 (tan 2 /? + 1 + A^ 2) ) 

- (c + 2 Cp ) - A A + im|(l - tan 2 0) + A« 

(23) 

where c, c M , ' , A A are radiative corrections, /i - 
Higgs mixing parameter, M A is the mass of CP-odd 
pseudo-scalar Higgs and mz is the mass of the Z boson. 

In fact these equations are just electroweak symmetry 
breaking (EWSB) conditions which are now solved for 
mi and 777,2. So from the above, we see that m\ and ra2 
can now be expressed in terms of \i and M^, which tells 
us that NUHM2 has the following free parameters: 

ra , rai/2, A , tan/3, /i, Ma, (24) 

where the trilinear soft breaking term Ao is not to be 
confused with the non-renormalizable term in inflation- 
ary scalar potential. 

In what follows, we want to find the regions of NUHM2 
which are compatible with the allowed mass range for the 
Higgs boson and the observed dark matter abundance. 
We will use two methods. One consists in identifying 
benchmark points which will satisfy all these require- 
ments, while the other method is more systematic and 
is based on a Markov Chain Monte Carlo (MCMC) scan 
of the NUHM2 parameter space. 

A. Identifying benchmark points for Neutralino 
dark matter 

To find interesting benchmark points, we use the 
micrOMEGAs code [29], coupled to the Softsusy spectrum 
calculator [30] and impose the following requirements: 

• The LSP must be a neutralino, 

• The relic density of the neutralino must be com- 
patible with the measured dark matter abundance 
by the WMAP experiment 0.1088 < tt D Mh 2 < 

0. 1158 pQ, 

• The LEP2 bound on the mass of chargino must be 
satisfied. It is given by > 103.5 GeV [54] . 

• The mass of the lightest Higgs must be within 
the range that is not excluded yet at the LHC, 

1. e. [115.5,127] GeV [63] 164] and more precisely 
equal to either m^ = 119 or m^ = 125 GeV. 

We scan the parameter space over the following range 
< ra < 3000 GeV and < m 1/2 < 2000 GeV (except 



tanp=10, A =-2m , |J=1000 GeV, M A =1000 GeV 




500 1000 1500 

m 1/2 [GeV] 



FIG. 3: (m ,m 1/2 ) plane for the NUHM2 model. We ex- 
plore specific configuration for tan/3 = 10, Aq — — 2mo, 
a = 1000 GeV, M A = 1000 GeV. Red strips are where 
0.1088 < QdmH 2 < 0.1158. Black lines show the two Higgs 
mass bounds. 



for cases where we require the Higgs mass to be about 
125 GeV Higgs as this pushes the upper bound on mi/ 2 
to mi/2 — 4000 GeV) and choose specific values of \i (the 
Higgs mixing parameter), tan/3 (the ratio of the VEVs 
for the up and down type fields) and Ma (the mass of 
the pseudoscalar Higgs). 

We set the mass of the top quark to the Tevatron value, 
i.e. m t = 173.2 GeV [50] and use the latest values of the 
following branching ratios BR(B S — ^ < 4.5 x 10 -9 

[66] and BR(b 57) = (3.55 ± 0.26) x 10~ 4 [49]. 

None of the scenarios that we find below can explain 
the measured value of the anomalous magnetic moment 
of the muon (g — 2) M ; the additional contributions in this 
model are indeed too small [24]. In what follows, we 
will assume that as long as the contribution of a given 
scenario is not greater than the measured value, the con- 
figuration is valid. 

The same observation and assumption are made when 
we consider the branching ratio BR(B+ — » r + ^ r ), know- 
ing the latest best average BR(B + — > t + v t ) = (1.67 ± 
0.39) x 10- 4 09] 

Our first results are summarised in Fig. In this 
figure the regions of the parameter space where the neu- 
tralino relic density is in agreement with the WMAP 
observations is represented by a red strip. The regions 
where the LSP is not a neutralino but a stau are coloured 
in green and the region excluded by the LEP2 limits on 
the chargino mass are represented in blue. The grey re- 
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gion corresponds to non physical configurations (in par- 
ticular, we find that the stop is tachyonic in most of this 
region) . 

Since we are looking for points which satisfy both the 
Higgs and dark matter constraints, we define benchmark 
scenarios as the points which lie at the intersection be- 
tween the red line representing the relic density and the 
two black lines corresponding respectively to a Higgs 
mass of rri} h = 119 and rri} h = 125 GeV. 

Our first conclusion for this choice of parameters is 
that it is hard to accommodate the correct LSP relic 
density with a Higgs mass of m^ = 119 GeV. There is 
only a small overlap when itiq = 248 GeV and mi/ 2 = 
834 GeV (denoted by 'a') corresponding to m % o = 351 
GeV and neutralino-stau co-annihilations. Indeed, be- 
low rai/2 ~ 830 GeV and for mo < 500 GeV, the stau 
and neutralino are almost mass degenerated. Hence the 
neutralino relic density mostly relies on neutralino-stau 
co- annihilations. For heavier neutralinos, both the coan- 
nihilation and annihilation rates decrease. As a result 
the LSP relic density becomes higher than the observed 
value. 

Unlike the case for m^ = 119 GeV, we find configu- 
rations intersecting the relic density and the rrih = 125 
GeV line. This assumes however that mo > 500 GeV 
and mx/2 > 800 GeV. In this region, the correct LSP 
relic density is achieved through CP-odd Higgs s-channel 
self-annihilations. To explain the observed abundance, 
the neutralino mass must be close to (but not exactly 
on) the resonance region. This leads to the relation 
m % o « M^/2 and thus implies that the neutralino mass 
is about m v o w 500 GeV for M A = 1000 GeV. This re- 
gion is actually referred to as the funnel region. Between 
the two red strips, the relic density falls below the ob- 
served dark matter abundance because the annihilation 
process becomes resonant and reduces the relic density 
too much. In total, we thus identify four benchmark 
points when m^ = 125 GeV. They are given by mo = 
1897, 2668, 1847, 2897 GeV with m 1/2 « 1100 GeV and 
correspond to the benchmark points 'b', V, 'd', V. 

The situation is a bit different in Fig. Q where Ma is 
larger. In panel (a) we find only one intersecting point 
for rrih = 119 GeV which is given by m 1 / 2 =411 GeV 
and mo = 837 GeV (benchmark point T). In this re- 
gion, the neutralino is mostly a bino. Hence we expect 
its mass to be roughly equal to Mi, which is related to 
mi/2 via Mi « 0.42mi/ 2 at any RG scale. Therefore, 
this point corresponds to Mi « m % o = 169 GeV and the 
relic density is achieved through neutralino-stop coanni- 
hilation. Indeed the upper grey region just above this 
point, corresponding to mo > 400 GeV, denotes to a 
tachyonic stop. The latter becomes the NLSP at the 
edge of this region, thus leading to an acceptable neu- 
tralino relic density. Such a point is nevertheless very 
likely to be excluded by recent LHC searches [71 , even 
though the precise limit on the stop mass depends on the 
model that is considered and Ref. [71] assumed a gauge 
mediated scenario. 



tanp=10, Ao=-2m , |J=500 GeV, M A =2000 GeV 




500 1000 1500 2000 



m 1/2 [GeV] 

(a) 

tanp=10, A =0 GeV, |J=1000 GeV, M A =4000 GeV 




500 1000 1500 2000 2500 3000 3500 4000 
m 1/2 [GeV] 

(b) 



FIG. 4: (m ,mi /2 ) plane for NUHM2: Panel (a) is for A = 
-2m , [i = 500 GeV, M A = 2000 GeV. Panel (b) is for fi = 
1000 GeV, A = GeV and M A = 4000 GeV. Both panels 
share tan/3 = 10. The red regions show 0.1088 < Q^DMh 2 < 
0.1158. 



For rrih = 125 GeV, we can define two more benchmark 
points. Both have a large scalar mass mo, namely mo = 
1715 GeV (benchmark point 'g') and mo = 2556 GeV 
(benchmark point 'h') respectively with m 1 / 2 = 1150 
GeV. Using the above relation between Mi and m^, we 
find that Mi is larger than /i, which implies an increas- 
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ing fraction of Higgsino component in the neutralino [51] . 
Such a composition favours annihilation channels such as 
X1X1 W+W~ , ZZ, Zh as well as neutralino- char gino 
co-annihilation and thus explains the vertical red strip in 
Fig. [4^a). Since the mass of Higgsino-like neutralinos is 
primary sensitive to the Higgs mass mixing parameter /i, 
we find that m % o w 460 GeV for both benchmark points, 
since they share the same \i and m^. 

In panel (b) of Fig. Q, we do not find any configu- 
ration compatible with both the 119 GeV Higgs and the 
observed dark matter abundance for the selected values 
of the tan/3, /i, Ma : A) parameters. However, we note 
that, within explored range of mo, any scalar masses at 
the GUT scale is compatible with = 125 GeV and 
the observed dark matter abundance. This leads to sev- 
eral possible benchmark points. The reason is that in 
this region the neutralino is very Higgsino-like. There- 
fore there is a large degeneracy between the two light- 
est neutralinos and the lightest chargino; hence the t- 
channel exchange process and the neutralino-neutralino 
and neutralino-chargino coannihilation mechanism con- 
tribute significantly to the LSP relic density 9 . As the 
gaugino contribution becomes negligible, m % o becomes 
mostly sensitive to the //-parameter. As a result the LSP 
is slightly heavier than 1 TeV. 

A summary of these benchmark points with their pre- 
dictions for important observables is given in Table ft]). 



B. A broader scan of the parameter space 

In the previous subsection we have identified a set of 
parameters for which the Higgs mass coincided with ei- 
ther rrih =119 or =125 GeV, and simultaneously lead 
to a dark matter relic density compatible with WMAP 
observations pQ. We now want to check whether the 
predictions associated with these benchmark points are 
generic or not. 

We thus perform a more general scan of the NUHM2 
parameter space. We now want to identify the regions 
of the parameter space which lead to a Higgs mass 
within [115.5, 127] GeV and a neutralino relic density 
within the WMAP measurements, namely QjjMh 2 G 
[0.1088,0.1158] using WMAP 7-year + BAO +# mean 
value pp. 

For this purpose, we use a MCMC coupled to the 
micrOMEGAs code along the lines described in [62 j. The 
total likelihood function is computed for each point cho- 
sen in the parameter space and is the product of the like- 
lihood functions associated with each observable. Since 



9 In the configurations of the parameter space leading to degen- 
erate two lightest neutralinos and lightest chargino it is possi- 
ble to embed the inflaton in a model where inflation is driven 
by the MSSM Higgses using the superpotential form: W D 
^ Hl 2^-3 55 . In this paper we will not discuss this model. 



we are not interested in characterising how statistically 
relevant the points that we found are but want instead to 
determine the full range of configurations that are possi- 
ble, we will not account for the number of occurrence of 
a given scenario. The drawback of such a method is that 
we cannot determine how likely a region of the parameter 
space is with respect to the other parts. The advantage 
is that very small (fine-tuned) configurations are kept in 
the analysis. 

For the Higgs mass and relic density, we define the 
likelihood as a function C\ which decays exponentially 

at the edges of the [x m i n ,x ma:c ] range, according to 



L>\ \X) X m { n ^ XfYiax 5 ®) ^ 2cr2 if X <C X m i n ^ 

( x ~ x m a x ) ^ 

= e if x > Xmax 

— 1 for X G [Xqjii-m %max\- (^^) 

with a a variance corresponding to the width of the 
[ x min-> x max] range and x the observable which corre- 
sponds in that case to either the Higgs mass or the LSP 
relic density. 

For all the other observables, we will use two types of 
likelihood. 

• For an observable with a preferred value \i and error 
a, we use a Gaussian distribution £2 • 

Q-/x) 2 

C 2 (x,ii,cf) = e 2^2 m (26) 

• For an observable with a lower or upper bound (set 
experimentally), we will take the function £3 with 
a positive or negative variance a : 

C 3 (x, /i, a) = - (27) 

1 + e 

We assume flat priors for all the parameters considered in 
this paper, and immediately reject configurations where 
at least one of the parameters fall outside the specified 
range. Points for which the calculation of the SUSY spec- 
trum fail (i.e. when there is no electroweak symmetry 
breaking or there is the presence of tachyonic particles) 
or the neutralino is not the LSP are also immediately 
rejected. At last, we do not implement the limits on 
sparticle masses from the LHC, since the squark masses 
that we consider are above the present limits. LEP lim- 
its on sleptons (and squarks) are nevertheless taken into 
account in micrOMEGAs. 

The known constraints that we impose from Particle 
Physics are summarised in Table. |H| and the range that 
we consider for the different parameters are given in Ta- 
ble ([in]). 

In Fig. ([5]), we see that most of the scenarios found by 
the MCMC involve TeV scale values of mo and m^, but 
no real feature emerges from the plot. 

Nevertheless, as illustrated in Fig. ([6|, there is a very 
strong correlation between the mass of the LSP and that 
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'b' 


(1897;1093) 


0.112 


125 


B 


3.113 


3.305 


6.07 


473 


A 


2249 


1955 




3 


'c' 


(2668;1085) 


0.111 


125 


B 


3.130 


3.241 


3.96 


473 


A 


2925 


2709 




3 


'd' 


(1847;1161) 


0.113 


125 


B 


3.108 


3.330 


6.05 


503 


A 


2249 


1914 




3 


'e' 


(2897;1152) 


0.112 


125 


B 


3.130 


3.252 


3.46 


503 


A 


3165 


2939 


4 


(a) >f> 


(837;411) 


0.109 


119 


B 


3.090 


1.731 


36.1 


169 


t 


952 


855 


4 


(a) 'g' 


(1715;1158) 


0.111 


125 


0.695+0.315 


3.076 


3.092 


7.82 


465 


X + '° 


2140 


1787 


4 


(a) 'h' 


(2556;1140) 


0.110 


125 


0.71B+0.295 


3.080 


2.921 


4.57 


462 


X + '° 


2850 


2603 



TABLE I: Benchmark points considered in this study and associated predictions for important observables. The Figures 
which they are associated to and the dominant mechanism (f, t coannihilations, x + '° exchange, A-pole) for the relic density 
calculations are specified in the last two columns of the table. The mass of the inflaton is at low scale. 



Constraint 


Value /Range 


Tolerance 


likelihood 


m h (GeV) [63, 64 


[115.5, 127] 


1 


£i(m h , 115.5, 127,1) 


Q x oh 2 pQ 


[0.1088, 0.1158] 


0.0035 


£i(Q x o/i 2 , 0.1088, 0.1158, 0.0035) 


Relaxing constraint on Qnh 2 


[0.01123, 0.1123] 


0.0035 


d(Q x oh 2 , 0.01123, 0.1123, 0.0035) 


BR(b — > sj) x 10 4 BUSH] 


3.55 


exp : 0.24, 0.09 
th : 0.23 


C 2 (10 4 BR(b^ s 7 ),3.55, 
V0.24 2 + 0.09 2 + 0.23 2 ) 


- 2) x 10 10 ESI 


28.7 


8 


A(10 10 (^-2),28.7,-8) 


BR{B S -> /x + /0 x 10 9 [SS] 


4.5 


0.045 


£ 3 (10 9 Br(B s /i + /i-),4.5, -0.045) 


Ap 


0.002 


0.0001 


£ 3 (Ap, 0.002, -0.0001) 


D (NUHM2\ rTTxl 

K B +^ T+i?T { SM ) [b5J 


2.219 


0.5 


£3(^+^+^,2.219,-0.5) 


Z -+ XiXi (MeV) 


1.7 


0.3 


£ 3 (^^XiXi,l-7,-0.3) 


a e+e-^ X X !3 


1 


0.01 


Cs(a e+e -^ x o x o 3 x Br( X 2,3 ^ Z X °i), 1, -0.01) 


xBr( X 2,3 -+ Z X °i) (pb) [67] 









TABLE II: Constraints imposed in the MCMC, from [70] unless noted otherwise. 



Parameter 


Range 


mo 


]0, 4] TeV 


mi/2 


]0, 4] TeV 


A 


[-6, 6] TeV 


tan /3 


[2, 60] 




]0, 3] TeV 


M A 


]0, 4] TeV 



TABLE III: Range chosen for the free parameters in the 
NUHM2 model. 



of the NLSP, suggesting that the neutralino relic den- 
sity either relies on the co-annihilation mechanism or a 
t-channel exchange of the NLSP (or both). The NLSP is 
found to be mostly a chargino, a neutralino and a stau 
as obtained for the benchmark points 'a'/g'/h'. The A- 
pole resonance corresponding to the benchmark points 
'b','cydye' requires however a certain amount of fine 
tuning (precisely because it requires m % o ~ Ma/2) and 
is not the most represented configuration found by the 



MCMC. 

The predominance of scenarios in which charginos are 
mass degenerated with neutralinos can be understood 
by inspecting Fig. For the configurations with 

Ao = —2 mo, the Higgs mass rrih tends to exceed the 
upper experimental bound unless one decreases the value 
of tan p 10 . . 

For such configurations, the sparticle masses are gener- 
ally too large for the sparticle-neutralino co-annihilation 
channels to reduce the relic density significantly and both 
the neutralino and chargino have a significant Higgsino 
fraction, Fig. [9j As a result, the possible channels to re- 
duce the neutralino relic density either involve CP-odd 
Higgs portal annihilations or neutralino-chargino mass 
degeneracies. 



10 Note that we didn't computed the amount of ElectroWeak fine- 
tuning in our NUHM2 scenarios. It was shown, for instance in 



|72| . that some NUHM2 benchmark points wherein Ao 
give non-negligible EW fine-tuning. 
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FIG. 5: Plot of the allowed parameter space in the (mo, mi/ 2 ) 
plane. We use the likelihood of the points as colour code. The 
darkest points have the highest likelihood. However they may 
not be statistically significant. 



FIG. 7: Higgs mass in the (Ao,tan/3) plane. Light Higgs 
can be found whatever the value of the trilinear coupling Aq, 
provided that tan /3 is small. 




m y0 [TeV] 



FIG. 6: Mass of the LSP vs the mass of the NLSP, depend- 
ing on the nature of the NLSP. Only points with dominant 
likelihood were considered. 



B component in xj 




tan [3 



FIG. 8: Bino fraction in the (Ao,tan/3) plane. 



The exchange of a pseudoscalar Higgs is actually sig- 
nificant when m % o ~ (as found for the benchmark 
points 'b','c','d','e') but neutralino-chargino coannihila- 
tion or chargino t-channel exchange are dominant when 
the Higgsino fraction is very large. In fact, among the 
configurations with a non-negligible Higgsino fraction, 
the larger the bino fraction, the more favoured the A-pole 
since small neutralino couplings to the Higgs can be com- 
pensated by having m % o closer to Ma- The distribution 
of points depending on their bino fraction is represented 
in the plane (Ao, tan (3) in Fig. ([8]). Clearly scenarios with 
Bino-like neutralinos are under represented, illustrating 
how fine-tune they are. 

Finally, we see from Fig. ^ that heavy neutralinos 
with a mass m x > 0.6 TeV have a large Higgsino fraction, 
thus suggesting even more dominant coannihilations with 
charginos (or annihilations through chargino exchange) 
and resonant annihilations via the pseudoscalar Higgs 
when the neutralino becomes fairly heavy. Interestingly 



though, for most these scenarios, the value of the /i pa- 
rameter varies between 500 GeV and 1.5 TeV but the val- 
ues which correspond to the highest likelihood are about 
H ~ 1 TeV, which is indeed consistent with a large Hig- 
gsino fraction. 

We can now investigate the distribution of points which 
satisfy the constraints on the Higgs mass and the dark 
matter relic density, Fig. (10). The points with high 



likelihood are 'smoothly' distributed within the observed 
relic density and Higgs mass range. 

Relaxing the constraint on the dark matter relic den- 
sity and allowing neutralinos to constitute only a fraction 
of the total dark matter energy density does not change 
the above features. The main effect in fact is to allow 
lower values of fi and reduce the mass degeneracy be- 
tween the LSP and NLSP. However these degeneracies 
are still present and larger values of the Higgs mass still 
give a higher likelihood. 
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Total likelihood 




FIG. 9: Neutralino composition. Top panel shows the bino 
content vs the neutralino mass while the bottom panel shows 
the Higgsino fraction. The colour coding corresponds to the 
likelihood of these points. 
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FIG. 10: Neutralino relic density vs the mass of the Higgs. 
The distribution of points show that any Higgs mass is asso- 
ciated to a high likelihood. 
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FIG. 11: (^cm^) plane for LLe and udd flat direction infla- 
tons respectively, where tan/3 = 10, Ao — — 2mo, \i — 1000 
GeV, and Ma — 1000 GeV. The yellow region corresponds to 
the Higgs mass of 119 GeV and the green region corresponds 
to that of the 125 GeV. Blue region is the same as in Fig. Q 
- excluded because of the mass bounds on chargino. Red and 
Black strips show where the dark matter abundance match 
within 0.1088 < tt D Mh 2 < 0.1158. 



IV. INDIRECT DETECTION OF THE 
INFLATON AT LHC 

In the previous section we have verified the validity 
of our benchmark points and could measure how fine- 
tune they are with respect to other configurations. In 



particular, we have seen that scenarios with large scalar 
masses mo require small values of tan j3 in order to not 
exceed the upper limit on the Higgs mass and lead to 
scenarios in which the neutralino has generally a non- 
negligible Higgsino fraction. 

We can now determine the inflation energy scale and 
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mass of the inflaton for these benchmark points. We 
will follow a similar approach as in Ref. [15] in order to 
estimate the inflaton mass which is compatible with the 
temperature anisotropy of the CMB data. 



A. Inflaton mass for benchmark points 



tanB=10, A =-2m , U=500 GeV, M A =2000 GeV 

JHM2 LLe case 



In Fig. (11) we have mapped the regions of Fig. 
onto (0o, rrifj)) plane. Panel (a) is for the LLe case and 
panel (b) is for the udd case. The yellow region is where 
the Higgs mass is equal to 119 GeV, and the green one 
is for rrih = 125 GeV. The blue area correspond to the 
LEP2 bound on the chargino mass. Vertical lines rep- 
resent the inflaton mass for which the Higgs mass inter- 
sects with the WMAP relic density measurements. We 
used the red color for rrih = 119 GeV, and black for 
rrih = 125 GeV. The grey shaded region shows where 
the NUHM2 inflation can explain the CMB observations. 
From these figures we see that if we want consistent de- 
scription of 119 GeV Higgs mass with the dark matter 
density measurements, inflation must happen in a range 
roughly given by: </> « (1.8 - 2.6) x lO^GeV for LLe 
and 0o ~ (2.1 - 2.8) x 10 14 GeV for udd inflaton can- 
didates, corresponding to a mass of ~ 300 GeV and 
m0 ~ 380 GeV respectively. 

In the case of rrih = 125 GeV Higgs, inflation should 
happen around O ~ (4.8-6.8) x 10 14 GeV for the V and 
'd' benchmark points, see Table. [l| yelding « 1900 
GeV for the LLe scenario and a slightly heavier udd can- 
didate. Another two possibilities correspond to the bech- 
mark points 'c' and V, see Table [l[ For V point we have 
inflation happening in a range of 0o ~ (5.7—8) x 10 14 GeV 
with a mass of the inflaton being around ~ 2700 GeV 
and similarly for 'd' we have range of 0o ~ (6 — 8.1) x 10 14 
GeV with « 2950 GeV. From cosmological point of 
view, the heavier Higgs boson is, the more of the parame- 
ter space for inflation which become compatible with the 
CMB observations we have. In general for the udd infla- 
ton, we have a larger running than in the LL widetildee 
case, essentially because of the running of #3. However, 
it is hard to appreciate this running visibly by comparing 
Figs. 11 (a) and (b), because of the large range of we 



have plotted. 

In Fig. 12 we are mapping the points of Fig. [4] (a). 
Here the red shaded region corresponds to the allowed 
relic density, whereas the yellow and green regions are the 
same as in the previous discussion. The Higgs mass of 
rrih = 119 GeV again implies the lower scale for inflation 
- i.e. 0o ~ (3 — 4.2) x 10 14 GeV, with an inflaton mass of 
around « 860 GeV. However for a Higgs mass of 125 
GeV, we find two energy scales related to different 
values: O ~ (4.2-6.2) x 10 14 GeV when m « 1780 GeV, 
and 0o ~ (5 - 7.5) x 10 14 GeV when ra « 2550 GeV. 
We obtain similar conclusions for Fig. ( 12 ) (b), where the 
udd direction gives a slightly higher mass for the inflaton. 
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FIG. 12: (^cm^) plane for LLe and udd flat direction infla- 
tons, where tan/3 — 10, Ao — — 2rao, the ix — 500 GeV, and 
Ma — 2000 GeV. The yellow region corresponds to the Higgs 
mass of 119 GeV and the green region corresponds to that of a 
125 GeV Higgs. The red and black strips show where the dark 
matter abundance fall within 0.1088 < QdmIi 2 < 0.1158. 



B. LHC predictions and Inflaton mass 

Our previous scans of the NUHM2 parameter space 
have selected neutralinos with a high Higgsino fraction 
when the neutralino mass falls within the 0.6 and 1.2 TeV 
range. It is now interesting to check the prediction for the 



13 



stop mass depending on the inflaton mass at TeV scale, 
see Fig. ( 13 ). We find that in both inflation scenarios, the 
inflaton mass is above 500 GeV and is associated with a 
very massive stop. For the udd combination, the lightest 
stop mass is constrained to be within > > ra^/3. 
Scenarios with the lightest stops (namely < 2 TeV) 
may offer a chance to probe the NUHM2 parameter space 
and thus a mean to determine the inflaton mass. 



Total likelihood 




guish between the udd and LLe scenarios. In addition, 
we find that staus in both scenarios can be lighter than 
1 TeV, thus offering another possible window for probing 
this model at LHC. Discovering a relatively light stau at 
LHC together with a specific stop mass would constrain 
the parameters of the model and thus provide a determi- 
nation of the inflaton mass. 



> 




FIG. 13: The lightest stop mass m t - versus the inflaton 
masses for udd and LLe, see Eq. (Uk. 



m- [TeV] 



FIG. 14: The correlation between stau mass, mf 15 and the 
lightest stop mass, . The color coding corresponds to the 

inflaton masses for udd and LLe. 



Such predictions have to be complemented by other ob- 
servables, such as the stau mass Fig. ( fl4] ). The prediction 
differs depending on whether the inflaton correspond to 
the udd or LLe inflation mechanism. For the LLe case, 
one finds that scenarios with 'light' inflaton (i.e. with a 
mass lower than 2 TeV) correspond to staus lighter than 
2 TeV and stops lighter than 2-3 TeV. More generally 
there is a correspondence between the inflaton and the 
stau masses, whatever the value of the stop mass. This 
correlation between the stau and the LLe inflaton mass 
can be understood because the inflaton is of leptonic ori- 
gin. Similarly, for the udd case, the inflaton mass is re- 
lated to the stop mass but there is no constraint on the 
stau. Although such a feature can be easily understood 
given the nature of the inflaton, using LHC observables 
and searches for sparticles could provide a way to distin- 



Specific observables such as B s — > and b — > 

57 are also interesting to consider. In particular, in 
Fig. (15), one can see that most of the scenarios which 



fall within the observed range of the b —> decay rate 
lead to a relatively large B s — > branching ratio, 

basically within 3 x 10 -9 and 4.5 x 10 -9 . Some scenarios 
are nevertheless excluded (i.e. with a contribution larger 
than 4.5 10 -9 ). This provides additional scope for de- 
tecting such scenarios at LHC since most scenarios are 
within the sensitivity of LHCb [66] . 

Finally, for completeness, we display the expected 
spin-independent elastic scattering cross section associ- 
ated with these scenarios in a Xenon-based experiment. 
We juxtapose on this plot the limit obtained by the 
XENON 100 experiment [57] which is extremely robust 
regarding the relative scintillation efficiency L e f f at this 
mass scale [59] (even though it may be affected by astro- 
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FIG. 15: The branching ratios of B s — > and b — >• are 

shown with the color coding corresponding to the likelihood. 
The shaded region shows points within b — »• 57 experimental 
and theoretical error bars. 



physical uncertainties, see [60j |6Tj and uncertainties on 
quark coefficients of the nucleon) , as well as the predicted 
limit for the XEN0N1T experiment. 

As one can see, most of the scenarios presented in this 
paper regarding NUHM2 are well below the present limit 
set by the XENON100 experiment and cannot be con- 
strained for the moment. However the projected sensi- 
tivity for XENON IT indicates that it may be possible to 
probe NUHM2 parameters in the forthcoming future if 
not already ruled out by the LHC. 



XENON100 (2011) 




XENONU (2017 forecast) 



the dark matter candidate), the Higgs mass constraints 
from LHC, and the constraints set on the inflationary 
potential to match the CMB constraints. 

We have considered two inflaton candidates (udd and 
LLe) for which the 'high' scale of inflation 0o is inti- 
mately tied up to the low scale physics at the LHC scale 
via the RGE, and which are compatible with the ampli- 
tude of the perturbations, Sh = 1.91 x 10 -5 and the 2a 
tilt in the power spectrum 0.934 < n s < 0.988 [1 . 

We used two methods. One consisted in finding bench- 
mark points and the other one in performing a more com- 
plete scan of the parameter space by using a MCMC code. 
Our main conclusion is that for most configurations the 
udd inflaton appears to be 'fairly light' but still heav- 
ier than 1 TeV while the LLe inflaton can be as light 
as 500 GeV. In both cases however it is possible to find 
configurations in which both the staus and the stops are 
potentially within the reach of the LHC, thus indicating 
that sparticle searches at LHC could actually provide a 
mean to constrain the inflaton mass for some subset of 
the NUHM2 parameter space. Such constraints would 
have to be cross correlated with the measurements of 
BR(£? S — >• /i + and BR(6 — >> 57) since all the scenarios 
found in this paper have predicted values for these two 
branching ratios very close to the present experimental 
limits. Finally LHC constraints or potential hints could 
be enhanced by the results of the forthcoming dark mat- 
ter direct detection experiments such as the XENON1T 
experiment. 



As can be seen from Figs. ( 11 ) and ( 12 ), hints of a TeV 
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scale inflaton together with the precise measurement of 
the Higgs mass would actually narrow down the scale of 
inflation. Combined with the Planck satellite measure- 
ments which is expected to constrain the range of the 
spectral tilt with a greater accuracy, one should actually 
be able to pin point both the scale of inflation <pQ and the 
corresponding mass at the scale of inflation, thus pro- 
viding a window on extremely high energy physics which 
also complements the current observations from the CMB 
radiation. 

To conclude, it is possible to embed inflation within 
MSSM. This interplay between inflation and dark mat- 
ter provides an exciting prospect where inflationary 
paradigm can be tested by the Planck, LHC, along with 
direct /indirect dark matter detection experiments. 



FIG. 16: The expected limit from XENON1T detector on 
the neutralino cross section (spin-independent) with respect 
to the neutralino mass. 



V. CONCLUSION AND PERSPECTIVES 

In this paper, we search for the regions of the NUHM2 
(a variant of the MSSM with non-universal Higgs masses) 
parameter space which are compatible with the observed 
dark matter abundance (assuming that the neutralino is 
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